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Abstract A succession and silviculture model (ZELIG.CBA) for broad-leaved Korean pine forest of Qiangbai Moutain 
Area was developed based on the framework of ZEIJG model and characteristics of Broad-leaved Korean pine forests of 
Changbai area. ZELIG.CBA model consists four sub-models: growth model simulating annual increment of individual tree in 
forest, regeneration model simulating annual establishment of different tree species, mortality model simulating annual age- 
related and stress-related morality; and silviculture model simulating the forest response to different silviculture scenario. The 
validation of the ZELIG.CBA showed that the accuracy of the model for the forest growth was more than 95%. The succes¬ 
sion from clear cutting site simulating showed that the ZELIG.CBA model was stable for long term simulation. And selective 
cutting experiment showed that the optimal scenario for broad-leaved Korean pine forests was removal volume 30% com¬ 
bining with 30a rotation. 
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Introduction 

Gaps in forests were often formed due to big trees fall 
down. Regeneration and growth in the gaps increase be¬ 
cause of the improvement of light condition. The gap will 
be filled with new regenerated trees. Watt named this 
forest dynamics as the forest growth cycle. Since D. Bot¬ 
kin developed the first Gap model (JABOWA model) 
based on the forest growth cycle theory, the research on 
gap model became one of the most active directions in 
ecology because of JABOWA model’s success. H. Shug- 
art and his students made great contributions to gap 
model’s development. Most Gap models came from 
FORET modbl (H. Shugart and West 1977). 

The Gap model of broad-leaved Korean pine forests in 
Changbai Mountains developed in the past focuses on the 
growth, succession and response to global climate change. 
Since the sustainable management of broad-leaved Korean 
pine forests js a very important issue for log production, 
the purpose of this study is to develope a model for simu¬ 
lating the dynamics of broad-leaved Korean pine forests 
and its response to different management scenario. 

Study area 

The broad-leaved Korean pine forests distribute mainly in 
Changbai Mountain and Xiaoing’an Mountains. Changbai 
mountain range stretches along the boundary between 


China and North Korea. The study area extends from 
127°06' to 128°16' E and 41°21' to 42°25' N. Elevations 
range from about 500 to 1100 m above the sea level .This 
the region is affected by monsoon climate and features a 
temperate continental mountainous climate, cold in long 
winter and short cool in summer. The annual mean tem¬ 
perature ranges from 2°C to 4°C. Average annual precipi¬ 
tation is 630-780 mm and mainly distributed in July and 
August. The accumulated temperature of a 10 °C is about 
1500 °C days, and frost-free period is 100-120 days. 

Soil type is the mountainous dark brown forest soil. The 
typical broad-leaved Korean pine forest, which consists of 
abundant species and complex structure, is multi-storied 
and uneven-aged. The dominant tree species are Piiuis 
koraiensis, Tilia amurensis, Fraximis madshurica, Quer- 
cns mongolica, Ulmiis Japonica and Acer mono. The cov¬ 
erage of shrub layer is 0.4 and herb layer 0.8. 

ZELIG.CBA model 

Model structure 

ZELIG.CBA model, a computer model for simulating the 
dynamics and silviculture of broad-leaved Korean pine 
forests in Mt. Qiangbai, is an individual-based “Gap 
model” which was derived from the ZELIG model of 
Smith and Utban (1991) and ultimately from the FORET 
model (Shugart and West, 1977) and JABOWA model of 
Botkin el al (1972a,b). ZELIG.CBA consists of four pri- 


1 The study was funded by the Opened Research Station of Changbai Mt. Forest Ecosystems, The Chinese Academy of Sciences 
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mary subroutines: regeneration, growth, mortality and 
silviculture. ZELIG.CBA modified regeneration sub¬ 
model of ZEL1G adapted to characteristics of the broad- 
leavtfd Korean pine forest; added algorithms of under¬ 
ground water affecting regeneration and growth of forest. 

Like its parent and sibling models, ZELIG.CBA simu¬ 
lates forest dynamics by accounting the establislimcnt, 
annual diameter growth, and mortality of each individual 
tree in the simulated forest. The basic approach in model¬ 
ing each of these demographic processes is to begin with a 
maximum potential behavior and subsequently modify this 
potential accounting to the status of the individual tree in 
the context the resources available on the modeled plot. 


^ * min (A/’) (2) 

Where f’s parameters are separately modify factor of tree 
growth environment such as light, temperature, soil mois¬ 
ture and soil fertility. 

Available Ugh! function 

species response to available light was described by the 
function as: 

f,(S h ) = c l (l-e) CliS "- ry) (3) 


Diameter growth function 

the optimal DBII increment function is the same as ZH- 


LIG 
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Where, c v is a scaling constant, c 2 determines the rate of 
change in growth relative to change in sunlight, and c, is 
the light compensation point (where net growth is 0). S h is 
the light available at height h, the light extinction through 
a forest canopy as a negative-exponential decay following 
the Beer-Lambert I aw, as demonstrated by Monsi and 
Sacki(1953): 


Where b h b 2 , b } are parameters of Richards function; G is 
tree growth parameter. Table 1 includes main specific 
species parameters used in ZELIG.CBA. 

The actual increment of DBII is as 


Sh = S(f exp|-& * LAI (/l)| (4) 

Where S„ is incident light, and lAI(h) is cumulative leaf 
area index (niVni 2 ) above height h. Here k is a constant 
describing light extinction through the canopy; generally, 
k takes on values in the order of 030-0.50 (Montcith 
1973), and in ZELIG.CBA is set to 0.40 . 


Table 1. The species-specific parameters In ZELIG.CBA 



species 

^ NUI 

/a 

/cm 

G 

//.« 

/m 

b 2 

b, 

(day’C) 

DD„„ 

(day’C) 

L 

D 

M 

F 

Firms karaiensis 

400 

160 

507 

31.84 

-0.304 

1.034 

600 

2850 

2 

4 

2 

3 

Larix o/i 

crisis 

200 

80 

700 

34.22 

-0.0269 

0.994 

600 

3250 

4 

2 

1 

3 

Abies ne\ 

ohrolepis 

250 

90 

225 

30.50 

-0.0297 

1.062 

680 

3050 

1 

2 

1 

3 

Picea ko 

raiensis 

300 

no 

389 

30.50 

-0.0286 

1.052 

620 

2850 

2 

4 

2 

2 

Belula platyphylla 

100 

80 

860 

28.86 

-0.0347 

0.920 

600 

3250 

5 

2 

t 

2 

Tilia amurensis 

300 

120 

482 

28.86 

-0.0347 

0.920 

600 

3250 

2 

3 

3 

1 

Fraxinus mandshurica 

250 

100 

596 

24.92 

-0.0378 

0.937 

850 

2650 

2 

4 

2 

1 

Juglans mandshurica 

250 

8(1 

403 

24.92 

-0.0378 

0.937 

850 

2650 

3 

3 

2 

1 

Populus davidiana 

80 

80 

980 

26.67 

-0.0329 

0.913 

650 

3250 

4 

2 

3 

2 

Quercus mongolica 

4CK) 

140 

385 

26.67 

-0.0371 

0.978 

700 

3050 

4 

1 

3 

3 

Acer mono. 

200 

70 

380 

19.11 

-0.0483 

0.939 

700 

3050 

2 

3 

2 

2 

Vlmus japonica 

250 

90 

650 

26.67 

-0.0329 

0.913 

700 

3050 

3 

4 

2 

2 


Notes: A m „, is maximum age, maximum DBM and maximum height of species respectively ; G is specific species growth parameter; /> 2 , />, is 


parameter of Richard function; DD„,„ and Dd max are specific species minimum and maximum growing degree days; /,, D, M and F are class of the 
species shade-tolerance, drought tolerance, wet tolerance and soil fertility intolerance respectively. 


Temperature 

Temperature effects on tree growth are modeled in terms 
of growing degree-days . Species response is computed as 
a parabolic function relating annual degree-days to the 
minimum and maximum values for a species (Botkin et al. 
1972; Shugart 1984.Solomon et al.1984): 




4(DD - DD mi mo)(DD - DD) 

(DD nnw(/) — DD oiiii(r)) 


( 5 ) 


Wliere, DD is a growing-season degree-day sum above 
5°C 
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Soil moisture 

The effects of soil moisture on tree growth is modeled in 
ZELIG.CBA in terms of a drought-day index (with the 
ground-water table is lower than 1 m below the surface in 
the growth season) and wet index (with the depth of un¬ 
derground water table less than 1 m in the growth sea¬ 
son) . 

In the case of depth of underground water table less 
than 1 meter in die growth season, the response of trees to 
the soil moisture is based on simple empirically function 
such as: 

L(SWI) -[max (’ V '-~ 5H ' -,0)r (6) 

wt 

Where: wt is species’ wet tolerance class, three classes 
were divided in ZELIG.CBA; SWI is soil wet index which 
relate to the depth of soil underground water table. 

When the depth of underground water more than 1 me¬ 
ter in the growth season, the response of trees to the soil 
moisture is a function as: 

f^DRI)-lmM( drl - DRI ,0)| OJ (7) 

drt 

Where: drt is species’ drought tolerance classes (the main 
species in Changbai Mountains are divided into five 
classes based on the characteristic of species); DRI is soil 
drought day index which was defined as the proportion of 
days during the growing season for which there is inade¬ 
quate soil mpisture (soil moisture below wilting point, -15 
bar). ZELIG.CBA calculate bodi 20cm depth soil drought 
day index and 100 cm depth soil drought day index. The 
20 cm depth soil drought day index is suitable to modify 
regeneration, and 100 cm depth soil drought day index 
suitable to modify trees growth. 

Soil fertility: 

The response of trees to soil fertility is expressed as expo¬ 
nential function such as: 

f S f ( SF ) - «i(l- exp[a 2 (SF -a 3 )]} (8 ) 

Where a b a v a 3 is empirical parameter which depends on 
the species; SF is relative soil fertility. 

Regeneration 

There are many micro-environmental factors involved in 
this process, It is difficult to simulate the mechanism of 
regeneration. Therefore we would simplified some proc¬ 
esses of regeneration. We assume in ZELIG.CBA that tire 
regeneration occurs if two conditions are satisfied. One is 
viable seed exist and the other is the seeds are dispersed to 
an area where environmental conditions are suitable for 
germination and seedling growth. The actual number of 


species regeneration every year is equal to the potential 
maximum number of species regeneration which is modi¬ 
fied by soil temperature and soil moisture. We express the 
process with the function as: 

N = * L * > /si, 1 } (9) 

Where f M andare soil moisture modifing factors above 
20 cm depth, /, is temperature factor, N m „ is species 
maximum regeneration number which relate with light in 
forest and species biological characteristic, ZELIG.CBA 
use the algorithm of KOP1DE (Shao,1994). 

Mortality 

Mortality is modeled as a stochastic event, and may arise 
from two sources: natural (age-related) mortality and that 
due to stress from site factors or suppression. The simula¬ 
tion of natural mortality related to aging is based on two 
assumptions: (1) 1% of individuals survive to reach a 
spccies-specific maximum age, and (2) mortality is con¬ 
stant with respect to age. These assumptions provide for 
suitably long life spans that the annual mortality rate can 
be approximated as: 

p .gr = 4-065 (10) 

Where 4.605 is the natural log 0.01 and A„,„ is maximum 
age in years. 

'Hie stressed individuals are those fail to achieve 10% 
of their potential growth increment or an absolute diame¬ 
ter increment less than 0.3mm for two or more successive 
years in ZELIG.CBA. It is assumpted that 1% stressed 
individuals survive to reach 5a for birch and poplar, 15a 
for Pinus koraicnsis and 10a for other species. Mortality 
by stress can be approximated as: 

4.065 M m „ (11) 

Where A„ re is the maximum surviving years of stressed 
individuals. 

Silviculture 

The silviculture subroutine consists of a series of regula¬ 
tions and algorithms, such as forest logging (selective 
cutting and clear cutting), tending and reforestation. Se¬ 
lective cutting algoritlim includes selective density regu¬ 
lation of choice cutting species, size and cutting order. 

Model validation 

Two means are used to test ZELIG.CBA. One is short¬ 
term simulation to test reasonability of ZELIG.CBA; the 
other is long term simulation to test efficiency of ZE- 
I IG.CBA. 

ZELIG.CBA is fit with the data of 100 pemianent plots, 
which were survived in 1980, 1986 and 1989. The results 
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(see Tab.2) show the accuracy of the ZELIG.CBA for 
simulating broad-leaved Korean pine forest growth. The 
mean error of ZELIG.CBA is less than 5% except for 
species of Picea koraiensis, Abies nephrolepis. The 
maximum error of ZELIG.CBA is less than 10%. 


Table 2. The simulation error of ZELIG.CBA 


Species 

Mean error % 

Maximum error % 

+/- 

Larix olgensis 

4.8 

8.2 

+ 

Abies nephrolepis 

6.1 

11.2 

- 

Picea koraiensis 

8.2 

12.7 

+ 

Betula platyphylla 

2.8 

6.5 

- 

Tilia amurensis 

3.5 

9.3 

- 

Fraxinus mandshurica 

4.6 

10.2 

+ 

Juglans mandslwrica 

5.0 

12.3 


Populus davidiana 

3.6 

8.5 

+ 

Quercus mongolica 

2.5 

7.6 

4- 

Acer mono 

3.7 

9.2 

+ 

Ulmus japonica 

4.2 

8.3 

- 


The long-term simulation test is the simulation of suc¬ 
cession from bare ground. The result shows that Populus 
davidiana takes dominance first and then Betula platy- 
phylla, the relative important value (RIV) of Populus da¬ 
vidiana is 40 and RIV of Betula platypliylla is 40 at 30a. 
The RIV of Populus davidiana and Betida platypliylla 


reach the maximum at 30~40a. The dominance is replaced 
by Tilia amurcnsis and Acer mono after 70-80a. RIV of 
those species reach climax at about 200a, and then decays 
after that. The dominant species is replaced by Korean 
pine ( Pinus koraiensis) after 300a. Pinus koraiensis domi¬ 
nates for a long time. There is no tendency that it is re¬ 
placed by other species. The RIV ratio of the species is 
relatively constant. In the broad-leaved Korean pine forest, 
this result shows the broad-leaved Korean pine forest is in 
the climax stage. 

The simulating result shows the regulation of the broad¬ 
leaved Korean pine forest succession. At the beginning of 
succession from bare ground, the first species invaded into 
the bare ground is the fast-growing shade-intolerant 
broadleaf species such as Populus davidiana and Betula 
platypliylla. After canopy closure, the relatively shade- 
tolerant species such as Tilia amurensis, Acer mono and 
Pinus koraiensis invade. When the shade-tolerant species 
growing up, the shade-intolerant species died out. Finally, 
Pinus Koraiensis dominate the forest for long time. For 
the long-term simulation by ZELIG.CBA, simulating re¬ 
sult is agreement with regulation of succession of broad- 
leaved Korean pine forest. Validation for long-term simu¬ 
lation shows that ZELIG.CBA is suitable for long term 
simulation. 



- Pinus koraiensis 

- - -Betula platyphylla 

- - - - Tilia amurensis 

- Populus davidiana 

- Acer mono 


Fig. 1. Forest dynamics after clear cutting 


Selective cutting experiment by ZELIG.CBA 
In the Experiment, 4 scenarios with rotation 30a compa- 
nied with removal 30%, 40% and 50% in wood volume 
were t*ken. Another scenario is contrast (undistuibed 
forest), simulating time is 300a. The results were showed 
in the Table 3 and Fig. 2 

The Wood volume of the experimental forest decreases 
temporarily after cutting and increase up to 250-300 


m 3 /hm 2 before next cutting when removal is 30% and 
40%in volume. There is increasing tendency in volume 
with removal 30% and 40% but undisturbed there is no 
increasing tendency with removal 50% in volume. And the 
remain volume in forest with 50% removal is only 90 
m 3 Jim' 2 . Therefore the scenario of removal 50% in volume 
is unsuitable for sustainability of broad leaved Korean 
pine forest and the scenario of removal 30% and 40% in 
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volume is suitable for sustainability of broad-leaved Kore¬ 


an pine forest. If we just consider wood volume. 



Fig.2 Wood volume dynamics with rotation 30 year 


Table 3. The DBH Distribution after Logging in 30a rotation for selecting cutting system 


Removal 

Volume 

Time /a 




DBH /cm 





<10 

10-20 

20-30 

30-40 

40-50 

50-60 

>60 


30 

350.0 

116.7 

133.3 

16.7 

0 

0 

0 


60 

2450.0 

66.7 

33.3 

66.7 

0 

0 

0 


90 

3166.7 

533.3 

0.0 

16.7 

33.3 

0 

0 


120 

2466.7 

916.7 

0.0 

0 

16.7 

0 

0 

30% 

150 

2683.3 

1083.3 

83.3 

0 

0 

0 

0 


180 

2716.7 

1066.7 

183.3 

0 

0 

0 

0 


210 

1833.3 

1200.0 

166.7 

16.7 

0 

0 

0 


240 

2266.7 

883.3 

200.0 

0 

0 

0 

0 


270 

2283.3 

1116.7 

166.7 

16.7 

0 

0 

0 


300 

2483.3 

1233 J 

250.0 

50 

0 

0 

0 


30 

216.7 

150.0 

100.0 

50 

0 

0 

0 


60 

2300.0 

100.0 

66.7 

16.7 

0 

0 

0 


90 

3883.3 

666.7 

16.7 

0 

0 

0 

0 


120 

2933.3 

1200.0 

0.0 

0 

0 

0 

0 

40% 

150 

1816.7 

1133.3 

50.0 

0 

0 

0 

0 


180 

1833.3 

1066.7 

66.7 

0 

0 

0 

0 


210 

2316.7 

916.7 

50.0 

0 

0 

0 

0 


240 

2800.0 

1033.3 

66.7 

0 

0 

0 

0 


270 

2466.7 

1466.7 

0.0 

16.7 

0 

0 

0 


300 

2166.7 

1733.3 

300.0 

0 

0 

0 

0 


30 

166.7 

150.0 

100 

33.3 

0 

0 

0 


60 

2066.7 

33.3 

0 

16.7 

0 

0 

0 


90 

3566.7 

433.3 

16.7 

0 

0 

0 

0 


120 

2566.7 

916.7 

0 

0 

0 

0 

0 

50% 

150 

2683.3 

683.3 

33.3 

0 

0 

0 

0 


180 

2433.3 

933.3 

33.3 

0 

0 

0 

0 


240 

2550.0 

1633.3 

166.7 

0 

0 

0 

0 


240 

2450.0 

866.7 

16.7 

0 

0 

0 

0 


270 

2883.3 

866.7 

16.7 

0 

0 

0 

0 


300 

2683.3 

1816.7 

166.7 

16.7 

0 

0 

0 
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The Table 3 shows the dynamics of DBH distribution in 
simulated forest and the minimum DBH for the culling. 
The minimum DBH for logging is about 30-40 cm and 
40~50 cm some about 20-30 an in the special logging 
year such as year 150,180 and year 240 when removal is 
30% in volume. The minimum DBII for logging is about 
20-30 cm some about 10-20 cm in year 120 and 240 
when removal is 40% in volume. Hie minimum DB11 for 
logging is about 20-30 cm some about 10-20 cm in log¬ 
ging year 120 and 210 when removal is 50% in volume. 
The optimal logging DBH for Finns koraicnsis and Picea 
koraicnsis are more than 30cm, and the optimal logging 
DBH of other deciduous species are 24-28 cm. (Wang and 
Xu 1980) In summery, the scenario of removal 30% in 
volume and rotation 30a is sustainable for Korean pine 
forest. 


Conclusions 

The validation of the ZELIG.CBA model shows that the 
accuracy of tlic model for the forest growth is higher than 
95%. For most tree species, the accuracy reaches 90%, 
that is reasonable for shot term simulation. Hie simulation 
for succession from clear cutting site shows that the ZE¬ 
LIG.CBA is stable for long term simulation. Therefore 
ZELIG.CBA is suitable for simulating growth and dy¬ 
namics of broad- leaved Korean pine forests in Changbai 
area. 

Selective cutting experiment shows the optimal scenario 
for broad-leaved Korean pine forests is removal volume 
30% combining with 30a rotation. But in the actual more 
than 40% wood volume has been removed combining with 
30a rotation in that area. Therefore, for sustainability of 
Broad-leaved Korean pine forests in ChangBai area, it 
should be reduced logging volume to 30% or the rotation 
should be extended. 
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